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Abstract 



Mathematical design of a quantum network with prescribed transport properties is equivalent to 
the Inverse Scattering Problem for Schrodinger equation on a composite domain consisting of quan- 
tum wells and finite or semi-infinite quantum wires attached to them. An alternative approach to 
study of transport properties of networks, based on solution of the direct Scattering problem, is 
also possible, but optimization of design requires scanning over the multi-dimensional space of es- 
sential physical and geometrical parameters of the network. We propose a quantitatively consistent 
description of one-particle transport in the network based on observation that the transmission of 
an electron across the wells, from one quantum wire to the other, happens due to the excitation of 
oscillatory modes in the well. This approach permits us to obtain an approximate formula for the 
transmission coefficients based on numerical data for eigenvalues and eigenfunctions of the discrete 
spectrum of some intermediate Schrodinger operator on the network in certain range of energy. 
The obtained approximate formula for the transmission coefficients permits to reduce the region 
of search in the space of parameters in course of optimization of the design. We interpret the 
corresponding approximate Scattering matrix as a Scattering matrix of the relevant solvable model 
which can be used when modelling the network by a quantum graph. 
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I. INTRODUCTION 



Modern interest to quantum networks is inspired by the engineering of quantum electronic 
devices and oriented to manufacturing of networks with prescribed transport properties. In 
simplest case of a single-particle, single-mode processes transport properties can be de- 
scribed in terms of Scattering. Quantum conductance was related to Scattering Processes 
in pioneering papers by Landauer and Biittiker, see [3, Q , and the role of the resonance 
scattering in mathematical design of quantum electronic devices was clearly understood by 
the beginning of the nineties, see [H 0, E^. Nevertheless practical design of devices up to 
now, see for instance 2^|, was based on the formal resonance of energy levels rather than on 
resonance properties of the corresponding wave functions. At the same time importance of 
interference in mathematical design of devices was noticed in j^, |27| and intensely studied 
3 0, 0, E3, see also recent papers 0, H H @ • 
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Experimental technique already permits to observe resonance effects caused by the shape 
of the resonance wave functions, see 0, El, EH- In our previous papers, see for instance 
OH 11 1, we proposed using the shape resonance effects as a tool for manipulation of the 
quantum current. In particular, we considered the process of resonance transmission of an 
electron in the switch across the quantum well from one quantum wire to another based 
on excitation of resonance modes inside the vertex domain- a phenomenon which was first 
mathematically described in [jjj ]. 

In actual paper we apply this basic idea to description of one-electron transmission in 
general quantum networks which can be formed on the surface of a semiconductor of quantum 
wells and quantum wires joining them or connecting them to infinity. The depth of the wires 
and the wells with respect to the potential on the complement of the network is assumed 
large enough to replace the matching boundary conditions on the boundary of the network by 
homogeneous Dirichlet conditions for electrons with energy close to the Fermi level E F = A 



in the wires, see The limitations we impose on the width 5 of the wires are practically 
not too strict. They are defined from comparison of the non-dimensional wave-number of 
electron in the wires on Fermi level A and the non-dimensional inverse spacing between 
eigenvalues of an intermediate operator on the resonance level A ~ A. The numerical 



analysis of the circular quantum switch, see ll| , shows that our approximation is applicable 



to the network constructed of a circular vertex domain with diameter d and quantum wired 
width 5 submitted to the condition d > 46. Hence our further assumptions on the width of 
the wires do not mean actually that the wires are "thin" or the connection of the wires to 
the quantum wells is weak, though we include below, lemma 3.1, also analysis of this case. 
We assume that the dynamics of the electrons in the wires is single-mode and ballistic on 
large intervals of the wires compared with the size of geometric details of the construction 
(the width of the wire or the size of the contacts). Our methods in this paper are based 
on analysis of the one-body scattering problem for the relevant Schrodinger equation on the 
network. 

A general Quantum network which is studied in this paper is a composite domain of a 
sophisticated form, a sort of a "fattened graph", see for instance [H^f. It is impossible to 
obtain an explicit solution of the Schrodinger equation on this domain in analytic form. On 
the other hand, the direct computation is not efficient for optimization of the construction of 



2 



the network with prescribed transport properties, because the resonance transmission effect 
is observed only at a multiple point in the space of geometrical and physical parameters. 

Substitution of the network by proper one-dimensional graph with special boundary con- 
ditions at the vertices looks like an attractive program aimed to obtainin g ex plicit formulae 
for solutions of the Schrodinger equation on the network, see Really, anal- 

ysis of the one- dimensional Schrodinger equation on the graph , or even on some hybrid 
domain as in [fj^, is a comparatively simple mathematical alternative, but estimation 
of the error appearing from the substitution of the network by a corresponding graph is 
difficult. Asymptotic behavior of discrete spectrum in composite domains with shrinking 
channels was studied in J43J see also references therein. Simultaneous shrinking joining chan- 
nels and vertex domains of a realistic compact network was studied in extended series of 
papers on the edge of last century. In particular in El, El the authors develop, based on 



[86j . a variational technique for discrete spectrum of Schrodinger operator on a "fattened 
graph" . They noticed that the (discrete) spectrum of the Laplacian on a system of finite 
length shrinking wave-guides, width e, attached to the shrinking vertex domain diameter 
e a , < a < 1 tends to the spectrum of the Laplacian on the corresponding one-dimensional 
graph but with different boundary conditions at vertices depending on speed of shrinking. 
In case of "small protrusion" 1/2 < a < 1 the spectrum in the corresponding composite 
domain tends to the spectrum on a graph with Kirchhoff boundary conditions at vertices; for 
the intermediate case a = 1/2 the boundary condition becomes energy-dependent, and for 
"large protrusion", < a < 1/2, the spectrum on a system of wave-guides, when shrinking, 
tends to the spectrum on the graph with zero boundary conditions at the points of contact 
so the vertices play the role of "black holes" according to jiif . 

In distinction from the papers 0, quoted above, our theoretical analysis of the 
Schrodinger equation on the "fattenTd ~aph» 

was based on a single-mode one-body 
scattering problem in the corresponding composite domain with several semi-infinite wires. 
The aim of actual paper is: an "approximate" description of transport properties of the 
network. We achieve this aim via derivation of an explicit formula for the Scattering Matrix 
in terms of the Dirichlet-to- Neumann map of a specific intermediate operator obtained via 
proper splitting ( see [sii ]) of the relevant Schrodinger operator £ on the whole network in 
a certain interval of energy. Construction of the system of intermediate operators {L A } is 
defined by the band structure of the absolutely continuous spectrum a a of the Schrodinger 
operator on the network the following way. 

The multiplicity n = n x of the absolutely continuous spectrum of the Schrodinger oper- 
ator on the network is finite on each finite interval. If the Fermi level E F = A is sitting on 
the spectral band A n = [T n , T n+ i] with the spectral multiplicity n = n(A), then the specific 
splitting £ — > £ A = L A Q)l A is applied to the Schrodinger operator on the whole network to 
create the Intermediate operator L M defined in the orthogonal complement of "open chan- 
nels" in the wires, see section 2. When matching the outgoing solutions of the corresponding 
intermediate Schrodinger equation with exponentials in the chopped-off open channels, we 
obtain scattered waves of the Schrodinger operator on the whole quantum network and cal- 
culate the corresponding Scattering matrix. The eigenvalues of the intermediate operators 
L A below the threshold T n+1 generically give rise to resonances of the Schrodinger operator 
on the whole network, see the formulae (f2"3J) and (|37)|) below, or become embedded eigen- 
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values. In simplest case,when just one quantum well with four equivalent wires attached is 
involved, see |ll| , and only one simple resonance eigenvalue A of the Intermediate operator 
is sitting on the first spectral band A 4 , close to the Fermi-level A, the Scattering matrix 
is presented approximately by the corresponding resonance factor (P). This factor contains 
the pole of the DAA-map at the resonance eigenvalue and the normal derivatives of the 
corresponding eigenfunction (p projected onto the corresponding four-dimensional entrance 
subspace E = E of the open channels,- the linear hull of the corresponding cross-section 
eigenfunctions in the wires: 

= _V^ + U< + 

1 ciThrp i \ ~ an, o 



— ^ — -+iK + (x;) 

■ en. >^^_^ +(A A ) 



(1) 



Here K + is the exponent obtained from oscillatory solutions in semi-infinite quantum wires 
and P+ 3 -^) (P+^ 1 is the resonance term associated with the resonance eigenvalue A 
and the corresponding eigenfunction (p . It can be presented in standard form || -P+^f || 2 Pq 
with the orthogonal projection P onto the one-dimensional subspace E°_ in E + spanned by 

p dtpp 
' dn 

The above "one-pole" approximation (0) of the Scattering matrix near the resonance 
eigenvalue A Q is obtained via substitution of the Dirichlet-to-Neumann map of the interme- 
diate operator by the resonance term only, neglecting all non-resonance terms and the con- 
tribution from the continuous spectrum. In case of the resonance triadic Quantum Switch, 
see , similar neglecting was possible in a small neighborhood of the resonance eigenvalue. 
The size of the neighborhood was defined by the "natural" small parameter « 1/20 which 
appears from above mentioned comparison of the inverse spacing of eigenvalues of the inter- 
mediate operator on the resonance level and the wave-number of the electron in the wire on 
the resonance energy. We derive also a "few-pole approximation" of the Scattering matrix 
and interpret it as a Scattering matrix of some solvable model of the network in certain 
range of energy. 

Actual paper has the following plan: after the formal description of the quantum network 
and introduction of an intermediate operator in the next section 2, we prove in section 3 that 
the intermediate operator L\ is self-adjoint, hence possess a resolvent kernel which has usual 
properties. Based on these properties we introduce the corresponding Dirichlet-to-Neumann 
map of the Intermediate operator and connect it with the Dirichlet-to-Neumann map of the 
Schrodinger operator on the Quantum well. The following section 4 contains the calcula- 
tion of the Scattering matrix and derivation of the corresponding "one-pole approximation" . 
Then in section 5 we suggest an interpretation of the "one-pole approximation" of the Scat- 
tering matrix on Quantum network as a Scattering matrix of a solvable model presented as 
a star-shaped graph with a single resonance node. In the last section we discuss briefly some 
historic details of the perturbation theory on the continuous spectrum. 

The aim of our paper is the description of an approach to the transport problem on the 
network based on the idea of the intermediate operator. This approach is demonstrated in 
details. But the total volume of this paper is not sufficient to discuss all subtle mathematical 
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details of each technical step. Hence in couple of places, which are clearly indicated, see for 
instance discussion of the singular spectrum in Theorem 3.1, we supply just a sketch of the 
corresponding proof. We plan to return to these questions soon in following publications. 



II. INTERMEDIATE OPERATORS 

Consider a network on a plane manufactured of several non-overlapping quantum wells (ver- 
tex domains) Q t , t = 1,2, ...T and a few straight finite or semi-infinite wires u m , I = 
1, 2, ... M of constant width 5 m and length d m < oo attached to the domains Q tl , ^t 2; • • • 
such that the bottom sections 7™, 7^ of the wire u m are parts of the piece-wise smooth 
boundaries <9f2 tl , dfl t2 , . . . of the domains Q tl , ^t 2 > • • • respectively. Each wire is attached 
to one (if the wire is semi-infinite ) or two domains (if it is finite) , so that the function m — > t 
or m — > ti,t% is defined although the inverse function may not be defined, since there may 
be a few wires u m connecting two given domains f2 tl , Vt t2 . In our further construction we 
use slightly extended vertex domains with small "cut-off's" 5uj m of the wires uT attached to 
them: SuT = {0<x<£: m },to avoid a discussion of the subtle question of Sobolev smooth- 
ness of restrictions of elements from the domain of the operator onto the bottom sections 
of channels joining the inner angles. Denote by 7 the orthogonal bottom sections of the 
"shortened" wires with cut-off's removed (from both ends, if the wire is finite) and attached 
to the corresponding vertex domains. We will use further former notations Q t , uj m for ex- 
tended vertex domains and for the shortened wires respectively. The previously described 
vertex domains and wires will not be used in the following text. We consider the spectral 
problem for the Schrodinger operator 

Cu = — A M u + V(x)u = Xu (2) 

on the network Q = Utf2 t U m u m with zero boundary condition on dQ. The kinetic term 
— A M containing the tensor [i of effective mass is defined as 



A ^ ) 1 j ;>- 



(2^ t ) _1 A, if e a 
if G uo n 



2u" 9x2 2a 1 - dy 2 

where fit plays the role of the average effective mass in the well Qt, x,y are the coordinates 
along and across the wire u m and 11 are positive numbers playing the roles of effective 
masses across and along the wire uo m respectively. We impose the Meixner condition at the 
inner corners of the boundary of the domain in form D(jC) C W% (Q). We assume that the 

potential V is constant on the wires V = V m and is a real bounded measurable func- 

tion V t on each (extended) vertex domain f2 t . Without loss of generality we may assume 
that the overlapping wires do not interact, otherwise we may treat the overlapping as an 
additional quantum well. The absolutely-continuous spectrum of the operator £ coincides 
with the absolutely- continuous spectrum of the restriction J2 m = m (t) °f ^ ne operator £ 
onto the sum of all semi-infinite wires u> m , m = m(t) with zero boundary conditions on 
the boundary \J m du m . This fact follows from the Glazman's splitting technique, ji^, see a 
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discussion in the proof of theorem 3.1 in the next section. The spectrum of the restriction I 
consists of a countable number of branches U 



lOO 



V m H — s ± , oo ) which correspond to the 



oscillating modes e™(y) exp±z y 2/z 
by the eigenfunctions of the cross-section e™(y) 



- 


S 2 TT 2 
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X 
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sny 



< A spanned 



2M <?m 



1,2, 



1,2, 



If the Fermi-level A does not coincide with any threshold, A ^ V m + 
then the linear hulls of eigenfunctions of the cross-sections corresponding to the open chan- 

2 2 

nels V m H — s ± - < A are considered as entrance subspaces of "open (lower) channels": 

m 2 ^m S rn 



E m (A) 
olds: 



^ B*<y/2vt(K-V m )6 m 
V s ^ >v /2^(A-y m )«5„ 



The entrance subspaces of "closed (upper) channels" 



E™(A), correspond to the upper thresh- 



+ V m > A. Denote by P™(A) the orthogonal projections onto the subspaces 

E±(A) C £2(7™) on the bottom sections of the finite and/or infinite wires uj m lying on the 
border of the (extended) vertex domain Q t - 

An important feature of the above quantum network is the "step-wise" structure of the 
continuous spectrum of the operator L with spectral bands of different spectral multiplicities 
separated by thresholds, see Jisj]. This property permits us to introduce the corresponding 
structure of intermediate operators. 

Note that elements from the domain of the Schrodinger operator C on the whole network 
Q belong locally (outside a small neighborhood of inner corners) to the proper Sobolev 
class D{£) C W^{QtJocl see 1471. For given Fermi level A > we define an intermediate 
operator, following 0, [ljj, EH where the suggested construction was introduced for star- 
shaped Quantum Switch. The intermediate operator is defined by the same differential 
expression as the Schrodinger operator £ defined on W -smooth functions in the orthogonal 
sum of spaces of all square- integrable functions on the components {uj m , Qt} of the network 
with special boundary conditions at the bottom sections 7" 1 G dQ t , t = 1, 2, ... of the wires 
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PT{A)u t 
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P™{A) 
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- Ut 






7 m 


7 m_ 



0. 



= 0, P™(A)u r 
?_ m (A) 







IT 



1 du m 


1 du t 




2f4n 971 


7t m 2/i t dn 


jr. 



0. 



(3) 



These boundary conditions are meaningful due to Embedding theorems, see |£1|, because 
functions from both W^i^lt, loc), W^iu}™) have boundary values in proper sense on the 
elements of the common boundary doo m r\dfl t - on the bottom sections 7 of the (shortened) 
wires. We also assume that the Meixner conditions are fulfilled at the inner corners for the 
component u t of the element in the vertex domains. The obtained split operator is denoted 
further in this section by £ A , and the whole domain of it in U t W / 2 2 (f2 4 , loc) © U m V4 / 2 2 



is 



denoted by D ( £ 



c^(n). 



. The Meixner condition can be presented as D y£ A 

To present the above boundary conditions in a more compact form we can introduce or- 
thogonal sums of upper (closed) and lower (open) entrance spaces ^ m © {E™) := E tt ± on 
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the wires attached to the domain We introduce also similar notations for the correspond- 
ing projections Pt,±u w of the vectors ^2 rn u m := obtained via restriction of functions de- 
fined on the wires attached to the given domain Q t onto the bottom sections 7™ G dVl t P\duj m . 
Similarly we can introduce the restrictions of the functions Ut defined in Q t onto the sum 
r* = U m 7 t m of all bottom sections 7™ of the wires u m attached to the domain fl t : 



u =u t 



and the corresponding projections Yl m P t ± = Pt ±1 an d then form an orthogonal sum: 



t,m 



Now the above boundary conditions can be re-written in a compact form as 

1 du n 1 du" 

2ji n dn 2jj} dn 



P + u n = P + u" = 0, P_ [u fl - u w ] = 0, P_ 



0. (4) 



Here fin, /J are diagonal tensors composed of values of the average effective masses in do- 
mains Q t and along the wires uj m respectively. The boundary condition 0] may be interpreted 

2 2 

as a "chopping-off" condition in all open (lower) channels 2 s ^± s2 + V m < A and a "partial 

2 2^ A 

matching condition" in all closed (upper) channels </ ±s. 2 + Kn > A. The restrictions Z m (A) 
of the split operator C onto the invariant subspaces TV? = V s 2 ff 2 , x L 2 (0, d m ) = 

— -j-^- + V / m <A 
2 2 

£7™ x L 2 (0, <i m ) which correspond to the (open) channels 2 S ± s2 + V m < A, are self-adjoint 

operators and admit spectral analysis in explicit form. We denote by l A the orthogonal sum 
of all operators I (A) in all open channels. The restriction L\ of £ onto the complementary 
invariant subspace is also a self-adjoint operator : 



A ' 



The above chopping-off construction can be applied either to open channels in all wires, both 
finite and infinite, resulting in the above operators L A and l A , or applied to the semi-infinite 
open channels only, resulting in corresponding operators L A , l A respectively, C A = L A © l A . 
One of the intermediate operators : La or L A can be more convenient for calculating the 
Scattering matrix, depending on the architecture of the network. In the remaining part of 
this paper we choose L A as an intermediate operator, chopping off only the open channels 
in semi-infinite wires. 

Some useful modification of the above chopping-off construction is also convenient: for 
given value A of the Fermi level one can choose the number A 1 > A and split the operator £ 
as L A © l A with l A , containing several closed channels, for given Fermi level, in the finite 
lower group. The formulae which are derived for the Scattering matrix in this case look 
slightly less elegant because of excessive number of matching conditions. But an essential 
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gain of this modification of the intermediate operator is the shift of the lover border of the 
continuous spectrum of the intermediate operator beyond the level A 1; at the minor price of 
appearance of a finite number of eigenvalues the intermediate operator L A below A 1 . This 
modification of the construction of the intermediate operator will be used in course of the 
proof of the Theorem 3.1 and further in section 5, see also a relevant "physical" motivation 
in the beginning of section 5. We postpone the standard derivation of the corresponding 
formulae for the Scattering matrix to the following publication, but concentrate further 
mainly on the contribution to the Scattering matrix from the discrete spectrum of L A . 



III. DIRICHLET-TO-NEUMANN MAP 

When studying Scattering on one-dimensional Quantum Networks one obtains the Scat- 
tered waves of the operator £ on the network via matching solutions of the split Schrodinger 
operator C A with elementary exponentials, or, generally, Jost functions, in the semi-infinite 
wires at infinity, see for instance 12, E^, 0, E3|- In case of more realistic two-dimensional 



wires matching of solutions from neighboring domains on common boundary is also an initial 
step of the perturbation procedure. 

Now we focus on the multi-dimensional techniques of matching based on a special version 
of the Dirichlet-to-Neumann map j^, |9(j modified for the intermediate operator on the 
Quantum Network. We need, first of all, the self-adjointness of the intermediate operator, 
the corresponding Green function and the Poisson map. Then, based on them, we suggest 
an algorithm for the construction of the Dirichlet-to-Neumann map of the intermediate 
operator which gives eventually the corresponding Scattering matrix. 

Most of facts collected in the following statement are very well known, but scattered 
in literature, beginning from 1841 Lot's of important analytic facts we need, and even in 
stronger form, can be found in |47L Esl]. 

Theorem III.l Consider the operators C A , L A , l A , C A = L A © I A , defined by the differ- 
ential expression the above boundary conditions 0), or just matching conditions on T, 
and Meixner condition at the inner corners of the boundary of the whole network. The 
domains of the operators consist of all properly smooth functions defined on the sum of 
all open semi-infinite channels TC + (A) = Ylm^+i-^) an< ^ ^ n ^he orthogonal complement 
7i-(A) = L 2 (f2) S m Wip(A) respectively. The absolutely- continuous spectra cr(/ A ), c(C A ) 
of the operators l x ,£ x coincide with the joining of all lower and upper branches in semi- 
infinite wires, respectively: 



a{l A ) = U 
a{L A ) = U 



2m± 4 



■<A-V„ 



S 2 7T 2 T/ 

+ V m , oo 



2 2 



>A-V„ 



+ V m , oo 



S 2 7l 2 



2/iM 

• m. 



The absolutely- continuous spectrum cr(C) of the operator C coincides with the absolutely- 
continuous spectrum of the the orthogonal sum C A = L A © l A : 

a(C)=a(C A ) = a(l A ) + a(L A ) 



8 



Besides the absolutely continuous spectrum, the operators C and L A may have a finite number 
of eigenvalues below the threshold of the absolutely- continuous spectrum and a countable 
sequence of embedded eigenvalues accumulating at infinity. The singular continuous spectrum 
of both C, C A is absent. 

The original operator £ on the whole network is obtained from the split operator C A = 
L A © l A via the finite- dimensional perturbation replacing the first of the boundary conditions 
HP by the corresponding partial matching boundary condition in the lower channels. 

Proof To verify the self-adjointness of the operator C A it is sufficient to prove that the 
operator C A is symmetric on the domain consisting of all properly smooth functions satisfying 
the above boundary conditions and it's adjoint is symmetric too. Then the operator L A is 
self-adjoint as a difference of self-adjoint operators £ Ql A . 

We introduce the notation U miS 7™ := T and consider the decomposition of the net- 
work Q = (U m uj m ) U (U s f2 s ). Then integrating by parts with functions continuously dif- 
ferentiable on each component of the above decomposition we may present the bound- 
ary form of the Schrodinger operator L without boundary conditions in terms of jumps 



u 



u 



u 



du 

On 



1 du u 



1 dun 



of the functions, and ones of the normal 



2fj} On 2/i t dn 

derivatives on bottom sections of the wires and the corresponding mean values {u} 



i 



2/i 



+ 



II dn 2u t dn 

(Cu, v) 
du 



dn 

du 

dn 



{v} " W 



(u, £v) -- 
<9v 



u 



dn 
dv 

dn 



dT- 



dT. 



Denote by P + the orthogonal sum of orthogonal projections onto the sum E of the entrance 
subspaces of the lower (open) channels in L 2 {T) and by the complementary projection 
, I = + P_. Inserting this decomposition of unity into the integration over T = U m)< 7™ 
we can see that the above boundary form is equal to 



p: 



p: 



p: 



du 

dn 

du 

dn 

du 

dn 



p A (du 

\ dn 



P + {v} " P+ {u} P, 

pt {v} - p k _ {u} p: 

P + A [v]-P + A [u]P + A 



P_ A fvl 



P A [ul P A 



d\ 

dn 

dv 

dn 

<9v 

dn 

on 



dT+ 



dT+ 



dT+ 



One can see from (jSJ) that the operator L supplied with the above boundary conditions 
is symmetric. Vice versa, if L + is the adjoint operator then the boundary form with u G 
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D(L), v G D(L + ) vanishes. Then P + {v} = 0, P+ [u] = follows from arbitrariness of values 
of the normal derivatives of the function -P+v on both sides of j™. Similarly another two 
boundary conditions J3J) for elements v G D(L + ) may be verified based on the arbitrariness 
of values of P_ {u} and P_ {§^}- Hence the adjoint operator is symmetric and thus self- 
adjoint. 

The restriction Z A of the operator C\ onto the lower channels is obviously a self-adjoint 
operator. We call it the trivial component of L A Hence the restriction L\ of it onto the 
orthogonal complement of lower channels- the nontrivial component of C A - is a self-adjoint 
operator too. The absolutely continuous spectrum a(L\) coincides with the absolutely- 
continuous spectrum of the restriction of the original operator £ on the upper channels 



mm ' 



7T 2 / 2 



Here is a sketch of the proof of this elementary statement. Consider a splitting of C with 
higher "formal Fermi level" A : >> A. The absolutely - continuous spectrum of the operator 
C is the same as an absolutely-continuous spectrum of the split operator, due to the finite- 
dimensionality of the perturbation: it is a sum <t(/ a ) of absolutely continuous lower branches 
corresponding to the trivial component and the absolutely continuous spectrum of L A . 

Hence it suffice to prove, that the spectrum of L A below the threshold A*. n is discrete. 

Consider inside the wires a solution of the homogeneous equation L A u = Xu with A < 
A < A r This equation admits separation of variables inside the wires 



1 d\ 1 d 2 u 



2/x" dx 2 2 f j, ± dy 2 



+ Vu = Xu, 



u 



i 



2/.i 



ETXI V 2M \| s' 2 2 M J 

c, sin ■ 



Then each locally smooth solution of this equation is exponentially decreasing in the wires 

V^T^IA A| 1,/2 

\u\(x,y) < C u e 1 .In particular, the Green- function of the operator L A is square- 

integrable, hence defines a bounded operator, if exists for given A < A : . Hence the spectrum 
of the operator L A does not have a purely continuous component on the interval (— oo, AJ. 
Assume that there exist an infinite sequence of eigenvalues {A r } which is convergent to 
A 2 < A. Without loss of generality one may assume that A r < + 2 1 for each r. The 
(normalized) eigenf unctions are continuous due to embedding theorems, 

sup \ cp r \ < CA, 

with some dimensional constant C, hence they are exponentially decreasing in the wires: 



\<P r \(x,y) <CAe " "" 1 
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due to maximum principle, and admit an uniform estimate of "the tail" in the wires 
\<p r \ {x,y)dm < C -j =k =. Then any linear combination <p = Yl c r l Pr °f the eigen- 

functions at least admits the estimate 

r 

due to the Parseval identity, and the estimation of the tail 

A 2 

\(p\ (x, y) dxdy < C 



( A l " A ) 

Now compactness of the unit ball in the subspace of all linear combinations of eigenfunctions 
ip r follows from embedding theorems on the compact part ft M = Q U jx 2 + y 2 < M j of the 
network and the uniform estimation of the tails. The compactness implies the fmiteness of 
the number of the eigenvalues of L A below A. 

The absolutely-continuous spectrum of the operator l A has constant multiplicity dinxE A 
on semi-infinite interval (A A ,00), where 

\ max 1 f ' 

^ T T 

A = max — — + V . 



m a -v 



21 ^ ^kP, 



2J1 1 - 5 



On the spectral band A. = ( A A , A A ) the spectral multiplicity of the absolutely contin- 

- 1 A \ max 1 min I x A " 

uous spectra of both operators £a, £ is equal to dimE^. In our case embedded eigenvalues 
are possible, but only a finite number of them on any finite sub-interval of the continuous 
spectrum, see for instance iE l3 • 

Singular spectrum of the Schrodinger operator with compactly supported potential on a 
finite joining of standard domains, is, generically, absent, see the discussion in 0], but the 
proof of this fact in special cases needs an accurate investigation. In our special case absence 
of the singular continuous spectrum may be derived from the fact that the partial boundary 
conditions imposed on the bottom sections of semi-infinite wires define the finite-dimensional 
perturbation of the operator C. Recently absence of singular spectrum in wave-guiding 
systems was discussed in EH by the method which may be applied in our situation too. 
Neverteless we supply below a sketch of possible proof based on classical ideas. 

Basic definition of if -smoothness in |82j |. volume 4, XXXIII. 7 admits a local formulation 
for the operator A, R X T~L C D(A) in Hilbert space H with respect to the operator H, R x = 
(H — XI), see 19711 . Conditions of local smoothness are actually derived in the theorem 
XXXIII. 25 of jS2|. They imply the corresponding local conditions of absence of the singular 
spectrum of H + A on the interval (a, b), if we know that operator H does not have singular 
spectrum on the interval and succeed to verify the local smoothness of the perturbation A. 
The deficiency elements used in course of the construction of the split operator for real A 
below the threshold A of the absolutely continuous spectrum of C satisfy the Helmholtz 

equation in the wires and hence decrease exponentially at infinity as exp [A A . n — A] 1/2 



x 
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Hence the finite-dimensional perturbation of jC a = L A © I A — > jC is relatively smooth below 
the threshold X min , which implies absence of the singular spectrum of H. 

The last statement of the theorem is obvious, see the corresponding calculations below 
in course of derivation the formula (}2"4*j) for the Scattering matrix in terms of the DN-map 
of the intermediate operator. This accomplishes the proof of the theorem. 



Note that accumulation of embedded eigenvalues of the operator £ to a finite point is 
discussed in j^, where also an example of the Schrodiner operator with exponentially 
decreasing electric field and accumulation of eigenvalues to the threshold is suggested. Note 
that presence of slowly- decreasing potentials in the wires implies, generally, an accumulation 
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of the embedded eigenvalues to thresholds. This general fact was discovered recently in 
In our case accumulation of eigenvalues on any finite interval of the spectral parameter is 
impossible because potentials on the semi-infinite wires are constant ( the electric field is 
compactly supported). 

Assume that the eigenfunctions of the absolutely-continuous spectrum of the operator 
L A and it's Green function G A (x, s) are already constructed: 

- A fl G + VG = XG + S(x-s). 

Then we can construct the Dirichlet-to- Neumann map (DN-map) of the operator L\. The 
standard DN-map is described in [fjil |9(| . The modified DN-map of the operator L\ can 

be obtained via projection onto the lower channels of the boundary current 2~-P+§^ 

of the solution of the homogeneous equation lu — Xu = with the boundary condition 

u = ur E E\ for regular values of the spectral variable A. The solution u tends to zero 
r 

at infinity if OA 7^ 0. It may be constructed using the resolvent kernel G A (x,x',X) of the 
operator L A . The corresponding Poisson map is given by the formula: 



Jt dn s 
Then the corresponding DN-map of the operator L A is 

1 dV A u 1 f 1 d 2 G A (x jS ,X) 



dG A (x, s, A) 



2n a dn 



f 1 d 2 G A (x,s,X) r A 



The normal current a ^ w " 7 belongs to W\ , £ > 0, hence the formal integral operator 
J r ^-^^^-ur(s)ds is unbounded in W^ 2 E (r). However, since C W^ 2 £ {Y) the pro- 
jection of the normal current onto is bounded from W^ 2 e (T) to W^ 2 £ (r), hence 

w V = -^V r ^7^" r(s)dseH/ "" e(r) - 

VAf A : W* /2 ~ £ (T) — > Wl /2 - £ (T). 

According to the previous theorem 3.1 there are only a finite number of eigenvalues A^ of the 
operator L A in any compact domain of the complex plane {A}. Then the following spectral 
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representation is valid for the corresponding DN-map: 



VAf\\) = £ 



^ + 9nj; ' \ + <9n r 

(2/i n ) 2 A - A- 



Y I i \ 2 + dn \ + ans dp, (5) 
^ A-p " 



where the normal derivatives are calculated on T and the integration is extended over all 
branches of the absolutely-continuous spectrum of the operator La emerging from the upper 
thresholds * , 2 + > ^. The tensor /z n of the effective mass is diagonal, as defined in 

the previous section. Based on above analysis one can prove that P^m e W 3 J 2 e and the 

r 

integral and the series © are convergent in proper sense due to the presence of the finite- 
dimensional projection P'^ onto the entrance subspacc E'^_ C W^ 2 ' ■ sec a similar reasoning 

in The whole finite-dimensional matrix-function T>J\f K (X) is bounded in W^ 2 ~ £ with 
ooles of first order at the eigenvalues A r of the operator L A and cuts along the upper branches 



7T Z F 

TjZ, s 2 



+ V m , oo) of the absolutely-continuous spectrum. 



Practical calculation of the DN-map DM (A) by the above formula (j5J) requires knowing 
of eigenvalues and eigenfunctions of discrete and absolutely-continuous spectrum of the 
intermediate operator L\. 

Another expression for the DN-map may be also useful. This expression for T>M is 
obtained in terms of matrix elements of the DN-map T>M of the compact part Qq of the 
network with all semi-infinite wires having zero boundary condition on their bottom sections. 

The spectrum of the Schrodinger operator L on the "extended" compact part fl of the 
network with zero boundary conditions on the bottom sections of the semi-infinite channels 
is discrete. The DN-map of L is an operator with the generalized kernel 

VN(x s A) = 1 d 2 G (x,s } \) 
(2/zJ 2 dn x dn s 

VN{\) : Wz /2 - £ (dtt ) — ► W 2 1/2 - £ (d£] )- 
We consider the restricted PAA-map framed by projections on E = L 2 (T): 

P E VM(X)P E . 

Denoting by <f> r (x) the orthogonal projections onto E of the boundary currents of 
eigenfunctions of the operator Lq we obtain the spectral representation of the framed DN- 
map of the operator L by the formal series 



/ VeSZT = \^ 1 (Prix)) {(p r (s) 

(2/iJ 2 A -A? • 4"(2A0 2 A -A? 



which is convergent in weak sense on smooth elements from E 
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To derive the formula connecting the DN-maps of the operator L and one of the in- 
termediate operator we present the framed DN-map of Lq as a matrix with respect to the 
orthogonal decomposition L 2 (T) = E + © E_ into entrance subspaces of the open and closed 
channels. Based on the observation E + £ W^ 2 £ , we see that the matrix elements T>Af + _(X) 
are operators mapping E± . . . into E± . . . respectively (for regular A) 

P.VNP. = ( K (VN) P + ^ m) P -):J W » ^ ) . (6, 

2 ]2 

Consider the basis {eY 1 } of all entrance vectors of closed channels, + V m > A, and 

introduce the diagonal matrix K A with elements kf 1 = —\^2/j} . 2^ a 2 + V m — A, on the 
entrance subspaces E™_ of closed channels „ \ l x2 + V > A. The corresponding operator 

K A _{\) := K_(X) has a bounded inverse (Kj)'* : W\ l2 ' e -> W 2 3/2 ~ e , £ > for A below the 

minimal upper threshold A A . 

Let u = {ut, u m } be a solution of the Schrodinger equation on the network with the 
boundary data u T £ E + on the sum V of bottom sections of the semi-infinite open channels 
and matching boundary conditions in the upper channels 

L 8uq 1 du u 



P+u n = u r , P A [u n - Uo,] = 0, P A 



2/i n <9n 2/i" <9n 







and standard matching boundary conditions on all bottom sections of the finite wires. De- 
note by u_ the projection of the solution onto the entrance subspace E_ in the upper 
channels in the semi-infinite wires. 

We calculated above the DN-map of the operator L\ as a projection onto E + of the 
boundary current 

VAf u T = t^P+^t- 
2fi a dn 

of the outgoing [i^ solution u of the homogeneous equation £ A u = Xu with the boundary 
condition u\ r = u r £ E, and the matching conditions in upper channels: 

P_ [u n - uj = 0, P^— — - ^— j = 0. 

Due to (P) his gives the following system of equations: 

P + (VAf) P + u T + P + (VAT) P + u_ = VAT A u r , 

P_ (VAT) P + u T + P (VAT) P_«_ = K_u_. 
Eliminating w_ we obtain the following statement: 

Theorem III. 2 The DN-map VAf A of the intermediate operator L\ is connected with the 
DN-map (VAf) of the operator Lq on the compact part Qq of the network Q by the formula 

VAT A = P + A (VAT) P A + - P A + {VAT) P A _ pA 7 - A _ - P A (VAf) P A + . (7) 
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The obtained formula for VJ\f A can be convenient when considering an independent shrinking 
JHf| of vertex domains and the quantum wires of the network, if the local geometry of the 
network permits it. The DN-map of the Schrodinger operator on the system of channels is 
contained in (JJJ) in explicit form, i.e. as diagonal matrices of square roots: the negative on 

the real axis matrix K = diag j — y^uJ" yj 2 J^- l s 2 K» ~ ^} anc ^ P urer y imaginary matrix 

with a positive imaginary part (on real axis) %K + = i diag j ^/2/i£ yfx — -^rp. Kn| • 

Matrix elements of the DN-map of the Schrodinger operator Lq in the domain Qq are rational 
functions of the spectral parameter with singularities at the eigenvalues of the operator L . 

Remark Note that the denominator {T>M)_ —K- := D(A) in the above formula (J7|) can 
be conveniently transformed due to the invertibility of K- = — |iT_| on the Fermi- level: 



D(A) = P'l (VAf) Pi -K_ = \K_\ 1/2 I + \K^ L " (VAf)_ \K. 



-1/2 , | 1-1/2 



\k\ 112 . 



where the operator in the bracket is is acting in W 2 ^ A 6 and may be extended onto L 2 by 
continuity. 

Our nearest aim is investigation of its structure and establishing conditions of invertibility 
of the denominator D(A). 

We can present the operator DM via Hilbert identity, see with the corresponding 
Poisson kernels P H for large positive H: 

VAf(X) = VAf(-H) - (A + H)V + H V H - (A + H) 2 V + H (L - XI)' 1 V H . (8) 

The first summand acts as an operator from W 2 ^ 2 ~ £ to W 2 2 ~ e . Hence, substituting (jHJ) into 
the preceding formula for the denominator we notice, due to E + C W 2 , and invertibility 
of K : W 2 ^ 2 6 — > W 2 ^ 2 £ , that the operator \K-\ 1/2 (VAf)_ \K-\ 1/2 is presented as 

iK_r i/2 p>AT(-/f)p_v_r i/2 - 

-(X + H)\K_\- 1/2 P A _V + _ H V_ H P A _\K^\- 1/2 - 
(A + H y\K^ 1/2 p A _P + H (L - XI)- 1 V H P A \KA l/2 - (9) 

The first term in this expression is a bounded operator in W 2 ^ 2 °, and may be ex- 
tended as a bounded operator onto E = L 2 (T). Others terms are compact in 
E. Moreover, one can prove, following similar reasoning in [69|, that the operator 

1/2 -f- — 1 1/2 

\K-\ P-P H \L — XI) P H P \K_\ is trace-class operator in E = L 2 (T). 

We will use the decomposition of the last term into the series of polar summands, and 
separate the resonance term ^l^ , with O = P_-j^, obtained as a projection of the bound- 
ary current of the resonance eigenfunction tp onto the entrance subspace E A _ of the closed 

channels. Then, oi 
the decomposition 



— 1/2 

channels. Then, omitting the factors \K\_ we may arrange the summands according to 



V xr_ 1 0o)(0o \- 1 <Ps){<Ps _ 1 0o)(0o r 

" - v - (2^y x - a° + 2^ (2^y a - ao • (2^y x - a° + ■ 
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where )C is the contribution to T>M from the non- resonance eigenvalues \° s 7^ Aq. The norm 
of the contribution K , as an operator from W'^ 2 e to W^ 2 e (for each e > 0), may be 
estimated by the spacing p(Ao) at the resonance level Ao ~ A. We explore the invertibility of 
the denominator D(A) in two special cases : for shrinking networks and for thin networks, 
see comments below. 



Lemma III.l // the wires are shrinking as y — > 4r and the quantum wells are shrinking 
as x — ► — , but the Fermi level is kept constant during the shrinking, then the norm of the 

operator K K, on functions of the variables x', y' is estimated in LziX) as 

|| K: 1 K \\< C (i-) . (10) 

Proof Notice that in course of shrinking the role of the resonance eigenvalue is being played 
by various eigenvalues of the Schrodinger operator on the quantum well. Generically the 
resonance eigenvalue is simple and the operator K /C is a bounded operator at the reso- 
nance value of the spectral parameter. Since the DN-map is homogeneous of degree —1 and 
K = K admits an estimate by from above, see (fT^j) below, the estimate of the non- 
resonance term is obvious if the distance of the corresponding eigenvalue A s from the Fermi 
level E\ remains strictly positive in course of shrinking. This gives the required estimate 
of the contribution from all non-resonance terms with eigenvalues which do not approach 
the Fermi level E\ in the course of the shrinking. The number of others non-resonance 
terms, which correspond to eigenvalues approaching the Fermi level in course of shrinking, 
is finite. It is sufficient to estimate the contribution from a single non-resonance term. Con- 
sider the term which corresponds to the eigenvalue closest to the resonance eigenvalue 

£ n 

Aq = t 1 ~ A. Notice first that the shrinking of the normalized eigenfunction ip 1 of the 

6 n 

dimensionless Schrodinger operator is described by the formula 



X 



x' \ -1 



[/„>.!* (£H 

Then the normal derivative of the shrinking eigenfunction is transformed as 

dcp[ _ dtp, 1 



dn dn e 2 



n 



If the Fermi-level A is kept constant in course of shrinking, then the closest to the Fermi- 
level resonance eigenvalue is shifted to the spectral point Aq ~ Ae 2 Q of the "dimensionless" 
operator, with A dimension [mass] \ Then the spacing p'(A) on the resonance level is 

p(E A s 2 ) 

calculated as — — 2 — -. The projection P E is presented via multiplication by the indicator of 

the bottom-section, and the norm of it is homogeneous first degree, hence proportional to 
5" . Combining all these facts we obtain the following estimate for the contribution to the 
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DN-map from on the resonance lever from a single non-resonance term 



p fM\ (p M sun I ^7 ^12 « 

' i ri„ / \ ; l ,h> su Pr I V v 9 ! f__ 



s. 



This result is in full agreement with the fact that the DN-map is homogeneous of order — 1. 

UJ 

Then we have II /C„ ||< CV. 

II II— ^ 

n 

One can easily obtain the estimate for K_ (A) at the Fermi-level: 



|| (K (A))" 1 || w max * /^V — ( 12 ) 

m y 7T 

if 2/x m (V") max m |i? A — << 7T. Summarizing the estimates (jllll2J) we obtain the an- 
nounced statement. 

Remark Introducing the positive operator \K_\ := — K_ one can derive a similar estimate 
in symmetrized form: 

\\\K_ \~ 1/2 )C \K_\- 1/2 \\<Co( £ - 

with some constant Co- 

The shrinking of the network, with the constant Fermi level A, can be applied to each 
quantum well and each quantum wire separately. Then the DN map on the joining of the 
wells is a direct sum of the DN-maps of quantum wells fit, t — 1,2, . . ., which implies: 



/C || < max 



sup r | V ¥s\ 2 £ 




|2 



\K_rX\K.r 2 ||<max,fe SUPr ^ ( ^ ) := C I - ) , (13) 



and 

ii i k r 1/2 r i k r 1/2 \\< m QY . / 

with an absolute constant Cq depending on the shape of the network. If the shrinking of 
details of the network (the finite and infinite wires and vertex domains) is independent, then 
max in the preceding formula is taken over all contacting pairs of details fl t , u m and the 
maximal value of the ratio appears in the right side. 

Definition We say that the network is thin on the Fermi-level A, if the condition || 
\K_\~ 1 K a \K_\~ ' \\< 1 is fulfilled for resonance values of the spectral parameter near to 
the Fermi level: A ~ A. 

This condition is obviously fulfilled for shrinking networks, if ^ << 1. But practically 
for given network and fixed A this condition can be also verified sometimes, based on direct 



calculations. In particular for the quantum switch [ll| based a circular quantum well radius 1 
with 4 quantum wires width 5 attached to it (centered at the points 0, 7r) the condition 
is fulfilled if 5 < 1/2, due to presence of some "natural" small parameter. 

The following statement (Theorem 3.3) reveals the structure of singularities in the above 
representation (JJJ) for thin (but non necessarily shrinking) networks. 
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To formulate and proof the statement we need more elaborated notations. Separate the 
resonance term in matrix elements of the DN-map of the operator L$ presented as a matrix 
with respect to the basis E±: 



P A _{VM)Pl 



1 0o)(0o 



_ \0 /o.. 



A - A; 



(2/iJ A-A j£ 

s^O 



1 




1 


0o)(0o 


(2Mo 


) 2 A-Ao 


(2Mo 


) 2 A-A° 


1 




1 


0o)(0o 


(2/i n 


) 2 A-Ao 


(2M n 


) 2 A-A° 


1 




1 





(2/0' 



A-Af 



A -A? 



(2/lnJ 



A-A° 



with /C_ + = (/C+_) , and 
P* (DAT) P K _ = 



0o ) (0o 



KM. 



0o)(0o 



(2// n ; 



A- A 



Here /C ++ , K ± , /C are the matrix elements of the contribution K to the DN-map from the 

non-resonance eigenvalues. The operators fC ++ , /C + _ are bounded in W^ 2 ~\Y), /C_ + = /C+_ 

for real A and K, acts from W^ 2 (T) into W^ 2 ~° (T) the same way as K_ does. Then 

the expression (J7J) may be written as 



T>J\f 



4> + o){4> + o 



0o)(<A 



(2/,J 2 A-Ag 







(2/,J 2 A-A° 
J 



i 4> )(4>o 



0o)(0o 



(2/0 A - A > 



(14) 



The positive operator — K_ = \K_\ can be estimated from below by the distance from A to 
the lowest upper threshold 



p (A) = min 



2f4 



1T 2 P 



(-K u, u) > p II u ir 



+ K, - A, 



Theorem III. 3 If the network is thin on the Fermi-level A, then the pole of the DN-map 
T>M at the simple resonance eigenvalue Aq of the operator L on the compact part Qq of 
the network (the singularity of the first addendum of Ffy) is compensated by the pole of the 
second addendum and disappears as a singularity of the whole function V>M so that the 
whole expression is regular at the point Aq. A new pole appears as a closest to A zero 

eigenvalue of the denominator P^ [T>M) Pi — K_ v = 0. 



18 



Proof If the network is thin, || \K_\ 1/2 IC (X )\K_ | 1/2 ||< 1, then the operator 



K -K 



k 



\k r /2 d 



is invertible: k 

the above product (fUj) can be found as a solution of the equation 



K I 1/2 /C \K I 1/2 | I if I 1/2 . Then the middle term of 



0o ) (0o 



(2/i f 



A- A 



+ /C_ 



if 



u = k /--A 0o ) (0o j ^ /)> 

where P = (2/i n ) (A — A°) + (0 O , & 1 0o )• Zeroes of the function D coincide with singu- 
larities of the middle term in the above formula (J7J) for the V>M and hence coincide with 
the eigenvalues of the intermediate operator L A . 

Substituting that expression into (|14j) we notice that all polar terms containing the factors 
(A — Aq) compensate each other so that the sum of them vanishes : 



VAf 



(2/0 A " A o 



A -Ar 



~U 1 



(2/i f 



k 



r 



D A - Ar 



1 /C + _fe <j>~) (<f 
1 <f) + ) (/C+_£f 1 f 



(2/i f 



A - An D 



1 /C+.Af 1 ) (/C+_Af 1 0" 



(2/io) 
1 + - /C+_fc" 



~> ( 



+ 



(2/i f 



D 

k' 1 0") 



+ 



+ . . . 



(15) 



where the dots stay for terms defining the regular summand of VAf in a small neighborhood 
of the resonance. 



Remark 1 One can see that on the first step of the approximation procedure we obtain 
the pole of T>M at the simple zero of the denominator T> with the same residue <f ) (0 as 
PA/", in full agreement with physical folklore. The residue, with a small correction obtained 
on the second step of the approximation procedure, is given by 

+ — JC + _k ) (0 + — K, + -k . 



19 



In particular this means that the portions -P+^f of resonance eigenf unctions of the interme- 
diate operator La in the entrance subspace can be found via the successive approximation 
procedure (on the second step) as: 

» + - /C + _ k 0~, 

if the network is thin. Similarly the shift of the resonance eigenvalue A° to the zero A° of 
the denominator D 

D(X D ) = 

can be estimated in first order of the approximation procedure as 

A =X -^j( Vo ,\K_\ *,„> + ..., 
where the dots stay for terms estimated by powers of 

■ -i/2 . .. 1-1/2 



j^-y(vvl^-l 1^—11^- 1 V? > 



Remark 2 Note that the exponent K_ can be estimated from below in terms of the 
wave-number of the electron in the wires at the resonance energy A o 



p (X°) = min \l 2u" \ I An — VL, 

' m. TCX 



TT I 



■+V m <A 



TT 2 / 2 



In case of thin networks the wave-number exceeds the contribution to the matrix element 
T>M__ of the DN-map of Lq from the neighboring non-resonance eigenvalues. 

The eigenvalues and the eigenfunctions of the discrete spectrum of the intermediate oper- 
ator L A can be found either by the minimizing of the corresponding Rayleigh ratio, or from 
the corresponding dispersion equation, 

Pt (VM) Pt -K_u = 

involving the DN-map of the Schrodinger operator L with zero boundary condition on the 
border of the compact part Q of the network. 

Theorem III. 4 The eigenvalues of the operator L A may be found as vector zeroes (A, u x ) 
of the dispersion equation 

D(\)v x =VAf-(\)v x -K_(\)v x = 0. 
In particular for values of the spectral parameter between the maximal lower threshold \ A 

± J max 

and minimal upper threshold X m _ n the dispersion equation takes the form 

,-1/2 , ,-1/2 

v=\K-\ VAf \K-\ v, (16) 
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-1/2 



VM— \K- 



with the bounded operator-functions \K- 
formed to an equation 

u-V(H, \)v = 0, 

with a trace class operator 



-1/2 



in L 2 (T). It may be trans- 



where 



(\K_\f 2 I + (\K_\y^(VAf_(-H)-(H + \)V:V H )(\K_\y^ (\K_\) 



V(H, A) := (H + A) 2 k 1/2 V + H n x V H k 1/2 , 



\K_\+ VM— (-H) - (H + \)V + H V h 



-1/2 



-1/2 



and TZ X , V x are respectively the resolvent and the Poisson map of the operator L on the 
compact part of the network with zero boundary condition. The corresponding scalar equation 
may be presented in the form 

det [I - V(H, A)] = 0. (17) 

Proof The projections u r G E_ of the eigen-function u, L A u = Am of the operator L A 

onto the cross-sections T of the open channels should fulfill the condition |^ = K*(\)u r . 

r 

On the other hand the restrictions of u onto the compact part of the network fulfills the 



corresponding homogeneous equation, hence 



du 
(in 



T>M (A). Matching both data with 



the boundary conditions (@J) we obtain the dispersion equation: 

K A _(X)u r = P K VH P A _ (A) u v . 

Due to the invertibility of K A (A) below A A we obtain the first statement (IToT) of the theorem. 
The second statement requires the iterated Hilbert identity for the DN-map, see |68| |: 

VM (A) = VM (-H) -(H+ X)V + H V H - (H + XfV + H R x V H . 

Both second and third terms of the sum in the right hand side are compact operators in 
L2<T), see the remark after the theorem 3.2, if the compact part Qq of the network has a piece- 
wise smooth boundary with Meixner boundary conditions at the inner corners. Moreover, 
the third term, after framing by factors K is an operator with a finite trace. Then the 
operator V(H, A) has a finite trace too and the dispersion equation may be presented in 
form 

det [/ - V(H, A)] = 0. 



IV. SCATTERING MATRIX 

In this section we rephrase some results of previous section in terms of Scattering Matrix 
and then derive and interpret the corresponding "one-pole approximation". 

Components of the Scattered waves of the Schrodinger operator C in the lower channels on 
semi-infinite wires u m are presented by linear combination of modes: oscillating exponentials 
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combined with eigenfunctions of cross-sections in open channels and decreasing exponentials 
combined with corresponding eigenfunctions of cross-sections in closed channels: 



Vm 



+ 



t 2 1 2 I v <A 



E 



e « b m,n L/ k e 



-A IE 



(18) 



where z/™ are the components of the incoming plane wave, and the finite matrix S = {S^n} 
is a Scattering matrix — the main object of our search. We present the above Scattering 
Ansatz in the following short form (|19|) introducing the notations: Ylmi e 
E + , {S l ^ n } = S, {s^ fc n } = s and the diagonal matrices in E+, E A : 



2^1 



A — V m — 



7T 2 / 2 



for &~ + V m < A and 



K 



2/Xm4 / V m + 



7T 2 l 2 



A 



for 



7T 2 fc 2 



+ V m > A, then the Scattering Ansatz in the wires is: 



+ e Sza + e 



SZA 



(19) 



To derive the formula for the Scattering matrix of the operator £ on the whole network 
Q we should substitute the above Ansatz into the boundary conditions of matching with 
components ip n of the Scattering Ansatz on the vertex domains 



0, P 4 



dip 
dn 



0. 



(20) 



where M |r = - V n ] |r and [g] | r „ /( 

of the above Ansatz in the wires are obtained by straight 



i av 



1 9 V> n 



T— r - The components ^} 

orward differentiation. However, 
the corresponding connection between the values of the components ^ip n , -^f" j of the Ansatz 
in the vertex domains Qt is given by the DN-map: 



VMip, 



2/z n <9n 
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Combining these connections and using the boundary conditions (|2U|). we obtain an equation 
for the undefined coefficients of the Ansatz p8jl . 

In this section we follow this standard program assuming that the network consists of 
"extended" vertex domains constituting the compact part Qq, see section 1, and several 
semi-infinite wires uj m , m = 1,2, .. . width 5 m attached orthogonally to Qq at the bottom 
sections 7 m , T = U m j m . The potential V(x) of the Schrodinger operator 

Lu = — A u + V(x)u 

takes constant values V m on the wires oj m and is a real bounded measurable function on 
vertex domains. Assume that the DN-map of the operator Iq on each vertex domain , and 
hence on fl is already constructed and presented (formally) by the spectral series in terms 
of the eigenfunctions (p r and eigenvalues A r corresponding to zero boundary conditions on 

on : 

2W(g,*0 = E x- xt" • (21) 

s r 

Here T>M is a direct sum of DN-maps of vertex domains. As a result of this steps, we obtain 
below the explicit formula for Scattering matrix in terms of matrix elements of T>M with 
respect to the orthogonal decomposition of the entrance subspace E = E + © E , see (|23|). 

On the other hand, choosing the Fermi level A we can define the intermediate operator La 
and construct the corresponding DN-map VAf as described in the previous section. This 
object is more sophisticate, because it already contains an information on closed channels, 
where the partial matching is already achieved. Really in second case we calculate Scattering 
solutions of the operator L matching limit values on real axis of the spectral parameter A 
of square-integrable solutions of homogeneous intermediate equations with boundary data 
from E + on T - with Scattering Ansatz in open challels only: 

ip u = e~ iK+x v + e iK+x Sv. (22) 

The boundary data of the solution of the intermediate equation are connected via the cor- 
responding (intermediate) DN-map: 



1 _ p diPo 



2/i" 



dn 



Assuming that v = Yl m i Vmie-T E + we obtain from the matching conditions the following 
equation: 

TW A [/ + £] = [-iK+ + iK + S] 

and the corresponding representation for the Scattering matrix, see below (|24|) . Summarizing 
two approaches to calculation of the Scattering matrix, described above, we obtain the 
following statement: 

Theorem IV. 1 The following two formulae are valid for the Scattering matrix on the Net- 
work: 
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a. The formula in terms of the standard DN-map T*N of the operator Lq on the compact 
part Qq of the network: 



S(\) 



(2W) ++ - (VAT). 



-K. 



(2W) , + iK , 



{VM) ++ - (VJ\f) 



+- (VAT) -K. 



(23) 



and 

b. The formula for Scattering matrix in terms of the DN-map T>M of the intermediate 
operator: 

In both formulae the denominator is the first factor of the product. 

Note that both above formulae are equivalent due to connection between the DN-maps 
of the intermediate operator and DN-map of Dirichlet problem on the compact part of the 
network, established in previous section. Though the formula (|24j) contains more sophisticate 
object T>M compared with T>N ', it may be more convenient for calculating resonances, 
see [HI, and for the description of transition processes, because the Scattering matrix is 
presented ClS db combination of bounded operators. 

One can notice that leading terms in the numerator and denominator of the above ex- 
pression (}2"4*j) for the Scattering matrix near the resonance eigenvalue A are : the polar 
term 

rim ( t\ > . r)/n ( q\ 

0o 0*0) (<^o( s ) 



g<gp (a) \ ^p^ dip (s) 



dn x 



dn s 



and iK + (\), which is estimated from below on the spectral band A^, 

1 



max ' min 



K (A) 



< 



2/4 



A- A 



A 

ma x 



1 
k 



(25) 



Assume that the remaining non-resonance part of DN-map 



p 9ip (x) , 
+ dn x 



9<P ( S ) 



A- A 



is subordinated to iK + in a small real neighborhood (A — e, A + e) := A £ of the resonance 



eigenvalue, that is : the operator K + 1/2 T>M Q K ^ 



1/2 



IKI 



- K is estimated as 
< d« 1, 



(26) 



on thin or shrinking network, with d = d(e). This, in particular, means that the spacing on 
the resonance level is large comparing with e. Note that, due to analyticity of both functions 
K + , DAT at the resonance level the above subordination condition (f2T)|) can be extended to 

C A 

some complex neighborhood of the resonance eigenvalue A Q 
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Form the unitary on real axis A combination of the leading terms in the numerator and 
denominator 



AIT I M 



iK, - 



(27) 



A-A„ 



Under certain conditions this combination may serve a "one-pole approximation " of the 
Scattering matrix near the resonance eigenvalue. Zeroes of S sit in upper half-plane and 
can be found from the algebraic equation A = A^ + i{K *(A) <f> , <f> ). If 



\(K +( f) ,<f) )\ <£ 15 and \(K' + (f) ,(f> )\ < q, < 1, 



(28) 



then a single zero of it sits in the above complex neighborhood U and it can be found by 
successive approximations procedure. 



Theorem IV. 2 If both conditions 



are fulfilled, with d small enough, then on the 



some real neighborhood U of the resonance eigenvalue the following estimate is true 



sup 



S(\)-S (\) 



2d 3 



< __r |k 1/2 (A A )| \K- l "(\~ 
1 - d 2 1 + V ° ;M + v °- 



-1/2 , A N 



(29) 



and i/iere exzst a smae zero of a Scattering matrix near to the above zero of a one-pole 
approximation. 



Proof. With notations introduced above the expression (|24|) can be presented on A e as 



A — A " 



ii", 



A-A* 



(30) 



K 



1/2 , ,A j-,-1/2 



with K — K ' PAT K 



K 



d. The formula ()27)1 can be presented in a similar form 



. 7 ^ 1/2 o ) (^ 1/2 ^o 

5 (A) = K 1,2 % + ^J^ ^K 1/2 := K- 1/2 SMk 1/2 



%L - 



Then 



K 



-1/2 



D_(X) 



A^Z 



S(X)-S (X) 
K 



D_(X) 



(31) 



J + - KIT 1 



(32) 
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In above real neighborhood the whole expression is small if K is small. Really, denoting 



— 1/2 

if) = K + <p and denoting by P^ the corresponding orthogonal projection, P^ = P + — P^, 
we have : 



D , 



D 1 — D 1 



1 



P +P — 

2(^;V o ,0 o )(a-a^ 



p 



D 1 — D 1 ) D 



(K- + ^ , <P Y + (A - A*) 2 ^ 

2(ir;> , O ) 

- (^;> o ,0o>-i(A-Aj) 



p, 



In particular on real neighborhood 





< i, 




< 1, 


(d; 1 -d: 1 )d, 


< 















and 



I-KD 



< - — -j, hence on real axis we have the estimation of the inner bracket 

l—d ' 



in (J32)) by and the estimation of the difference (|32j) on the real neighborhood A e 

d 



S Sr. 



< 



1 - d 



sup |< /2 (A)| sup |ir; 1/2 (A)|. 



(33) 



A€A„ 



ASA F 



For small q and hence small e, the above estimate can be replaced by the 



S S n 



< 



5i< 2 (a:)i liC'TOi, 



1/2 - A ^ 



1 - d 



with some coefficient -B greater then 1, for instance B = 3/2. Similar estimation may be 

Q 

obtained inside the corresponding complex neighborhood U e on the small circle C e centered 
at the single zero A = A A + i(K + (\ )(fi , O ) of S , with radius equal (K + (\ )(f) Q , (j) Q ) . Then 
we have similar estimate, 



S S n 



< 



2d 
1 - d 



sup |i^ /2 (A)| sup \K '''(A'J 



1/2 / . A\ 



AGC F 



ASCV 



b\k':'(K)\K"(K)\ 



1 - d 



(34) 



where we replaced maximal values of norms on the circle by values of norms at the resonance 
eigenvalue A Q , assuming that e is small. If the condition 



2d 
1 - d 



s|< /2 (a a )| \k;-(\;)\ < i 



1/2 . A, 



is fulfilled, then there exist, due to the operator version of Rouchet theorem, [36], an isolated 
zero of the Scattering matrix inside the circle C e . 

Summarizing results obtained in actual and previous sections we see that the step-wise 
structure of continuous spectrum on the quantum network permits us to reduce calculation 
of the Scattering matrix on the spectral band A. = A A , A A to matching solutions in 

• A max 1 min '- 
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neighboring domains via proper Dirichlet-to-Neumann map. Hence the large-time asymp- 
totic of corresponding non-stationary processes is defined by the shape of eigenfunctions of 
discrete spectrum of the intermediate operator. One can use the variational method for the 
construction both the eigenfunctions and eigenvalues of the intermediate operator below the 
threshold A , and thus obtain the Scattering matrix from the solution of the variational 

min 1 J 

problem. 

The derived estimation of the deviation of the one-pole approximation from the Scatter- 
ing matrix and the estimation of positions of corresponding resonances are too crude. The 
following "physical" motivation shows, that possibly contribution from some non-resonance 
eigenvalues and from the continuous spectrum in the previous expression (j23| for the Scat- 
tering matrix may be neglected under much wider conditions. This motivation shows that 
the one-pole, or possibly few-pole approximations, may give a reasonably complete physical 
picture of dynamical processes on the network. 

Assume that A is the Fermi level in the wires and T is properly scaled temperature. Then 
one may obtain an approximate expression for the Scattering matrix in an essential interval 
(A — T, A + T) of energy by substituting into properly reduced expression for the DN- 
map ZW containing singularities only in that interval. Really, one usually estimate the 
number of excited modes in Fermi-systems with discrete spectrum , considering the corre- 
sponding essential interval, because dynamics of electrons outside the interval is suppressed 
due to the Fermi distribution. We just extend this principle to the case when continuous 
spectrum is present, and implement it via replacing of the DN-map of the intermediate 
operator by the "essential part" of it : 

p dip r (x) \ lp dipr(s) 

.nA^/W 1 + dn x I \ r + dn r . 



VM — > DAT (A) 



T 

"n 



(2/iJ 2 ^ . A-AA 



[ A-T <A A <A+T 



J2 T^-Ti / + _ + 9ns dp. (35) 

1/ I * 2 ' 2 - A ' ^ J(A-T< P <A+T) X P 

2fi 7n o m 

Substituting the the essential part of the DN-map into the expression (j23| we obtain the 
corresponding "dynamical approximation" for the Scattering matrix on the essential interval 
of energy: 

S(X) W -VJ^—lKl- (36) 

The reduced essential part T>J\f^ of the Dirichlet-to-Neumann map of the intermediate oper- 
ator may contain one or several polar terms and an integral on an interval of the continuous 
spectrum. We will consider in further text the case when the essential interval does not over- 
lap with continuous spectrum of the intermediate operator A + T < \ min and neglect the 
integral on the continuous spectrum. We call the essential interval of energy [A — T, A + T] 
an essential spectral band and denote it by A . One may assume, that for narrow wires the 
continuous spectrum does not overlap with essential spectral band. The resulting essential 
expression (j3fi|) for the Scattering matrix is actually a Scattering matrix of some solvable 



27 



model, similar to one constructed in 0, 3^, z2| • In particular we again obtain the "one-pole 
approximation" Q for the Scattering matrix near the resonance 



S(X) 



P7V A + iK A 
£W A — iK A 



^ dn x 1 x ^ on s 



+ iK+{\) 



~ on x 1 > ~ qn s 



-iff+(A) 



if only one resonance eigenvalue A Q ~ A of the the intermediate operator is sitting on the 
essential spectral band A . 



V. STAR-SHAPE MODEL OF A QUANTUM NETWORK 

The aim of this section is: to construct a solvable model for general Quantum network 
with a finite number of semi-infinite wires. The model network will be constructed of a 
single vertex, where the essential part of the intermediate operator is substituted by some 
finite-dimensional operator, and several quasi-one-dimensional quantum wires attached to 
it, emulating the open channels. We will choose the parameters of the model such that 
the corresponding Scattering matrix will coincide with the essential Scattering matrix of 
the Quantum Network (}3l)j) . Similar construction for the simplest Quantum Network - a 
three-terminal Quantum Switch - was developed in preprint [fij], see Fig.l below. We will 
use in this section notations similar to ones introduced in section 2 



but we assume now that vector functions U , U present elements of functional spaces not 
on the bottom sections of channels, but on the vertex domains and wires respectively. With 
these notations we can present the basic Schrodinger equation on the vertex domains and 
on the wires in symbolic form 



1 



AU + V U a = EU 



1 d 2 u" 1 d 2 u u 



, +V U U"=EU", (37) 
/i" ox ay 

with a bounded measurable potentials V n = {V t } on the wells and constant potentials 
V" = {V m } on the wires. Consider the the entrance eigen- vectors e of cross-section and 
corresponding orthogonal projections p™ in spaces of square-integrable functions on bottom 
sections of the wires u™ , see section 2, and form the entrance subspaces E ± for open and 
closed channels with corresponding projections 

P± = Y1 ® p 7' dim ^+ = N 
± 

with summation over the thresholds below and above the Fermi level A respectively. 
For instance: we take for E + the sum of all p with indices satisfying the condition 
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1 V 

2 — ^ — \-V m < K for E . For values of the spectral parameter on the minimal spectral band 



2m± S 2 

, \ L containing the Fermi level A, we introduce an operator-valued Wave-number 

%n 1 max ^ ' A 

for open channels positive square root from 



1 2 ; 2 
1 7T / 



{+} 



with summation {+} over all open channels, and the positive decrement P__ for closed 
channels: 



2 ,2 



p =£ kr— + P, =V -A. 

{-} k m m J 

The introduced operators are connected with operators K ± from the previous section 



K ± = V2//P ± . 

Corresponding exponential solutions of the Schrodinger equation in the wires are: 
e e, with e £ E + , and e , e with e £ E_. 

Depending on architecture of the Quantum Network one may construct the intermediate 
operator either blocking open channels in all wires, as suggested at the end of section 2, or 
blocking open channels only in semi-infinite wires. When constructing the solvable model of 
the network we also have a choice: either we present the whole network as a single star-shaped 
graph with a single resonance vertex, or as a joining of star-shaped graphs with resonance 
vertices, connected by finite quantum wires. In second case the parameters of the finite 
wires enter the model as explicit parameters. We guess that the second representation is 
more convenient for detailed engineering analysis of networks and especially for optimization 
of their transport properties. Nevertheless we choose now the former construction, firstly: 
as more elementary and, secondly: just to include in our construction non-trivial ( bent 
of periodic) internal finite wires. When analysis of the first construction is done, we can, 
at least in principle, develop analysis for the second one, with straight connecting wires, in 
algebraic way, based on star-shaped elements. Note that importance of star-shaped elements 
was noticed long ago and remains motivation of numerous recent papers, see for instance 
|2fii [3Qj | . We prove here that the star- shaped graph with a "resonance vertex" , see our 
construction below, is an universal object in spectral theory of Quantum Graphs . 

When chopping off the open channels in semi-infinite wires only, we obtain a self-adjoint 
intermediate operator L A and construct the corresponding DN-map VM and the approx- 
imation of the Scattering matrix on the essential spectral band, (jHfijl . The corresponding 
Scattering Ansatz (|2*2*|) in semi-infinite open channels satisfies the equation 

1 dU 

-71?- + vu =xu - (38) 
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We choose the outer ( "trivial" ) the part of the model - in open channels - in form of original 
matrix Schrodinger operator 

1 dU 

I U = — t-t^ + v u 

H ax 

in L 2 (E, [0, oo)) and we will use the same Scattering Ansatz (j22j) in open wires as for the 
original problem. But the construction of the vertex part of the model will be done with a 
major change of the original intermediate operator. 

Based on the above physical motivation we assume that only a finite number A T of 
eigenvalues of the intermediate operator sits on the essential spectral band and no change 
of multiplicity of the continuous spectrum occurs. Then we will substitute the intermediate 
operator on the Network by some finite-dimensional positive hermitian matrix A acting in 
an abstract space E A , dimi^ = N T . The eigenvalues and normalized eigenvectors of the 
operator A will be denoted by k s , e s correspondingly, s = 1, 2, . . . N T . The eigenvalues and 
the boundary parameters of the model j3, see below (J3Hj) . will be defined later based on 
comparison of the Scattering matrix of the model with the essential Scattering matrix. The 
resulting operator we obtain is a simplest star-shaped model of the Quantum network 
with the compact part substituted by a sophisticate Quantum Dot. In our construction we 
follow the receipt of construction of zero-range potentials with inner structure described in 
E3| , see also references therein. 




FIG. 1: Three-terminal Quantum Switch 



The corresponding model for the three-terminal Quantum Switch is used in [10] for estima- 
tion of the parameter regime of the switch. 

Restrict both / and A to symmetric operators onto the corresponding domains: I — > l — 

defined on functions vanishing near near x = 0; the operator is restricted A — > A — 



I 
A 



such that given subspace N 6 E A would play the role of deficiency subspace at the 



spectral point i : N = N dim N = d. Then dual deficiency subspace is N_i = ^^jA. 
Define the domain of the restricted operator A as = -j^j [E A AJ of A . It is not 
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dense, since A is bounded. Nevertheless, under condition that the deficiency subspaces 
N ± = N ±i do not have common elements, see for instance , the extension procedure 

for the orthogonal sum l © A can be developed with use of the symplectic formalism, see 
for instance [72|. The case dimE A = 1 is not formally covered by the above condition, 
but the corresponding extension also was constructed in 00 and the relevant fonnnlae 
for the scattering matrix and scattered waves remain true. The extension of the restricted 
operator A is reduced to selection of a Lagrangian plane of the sum of the corresponding 
boundary forms, see for instance [H3, Zil, and below, see E21 where the Lagrangian plane 



is chose based on proper boundary condition. The same description of extensions remains 
valid in case dimE A = 1. We will use this fact below, see EH 

The boundary form of the differential operator is obtained via standard integration by 
parts: 

J t (U, V) = (l + Q U, V) - (U, l + Q V) = (U'(0), 1/(0)) - (17(0), ^(0)), (39) 

where £7(0), V(0) G E + and the derivatives are taken in positive direction. 

We assume that the Hermitian operator A : E A — > E A , dim E A = N T , is defined by its 
spectral decomposition (on discrete spectrum) 

s 

2 

where k s are positive eigenvalues and P s are the corresponding spectral projections, P s = 
e )(e a . Choose an ortho-normal basis in Ni {ft} = {ft 7 i} , I — 1, 2, . . . , d, d — dim N < N T , 
as a set of deficiency vectors of the restricted operator Aq. Then vectors jfzpjfi form an 
ortho-normal basis in the dual deficiency subspace iV_ r We assume that the subspaces Ni 
and N_i do not overlap, hence their (direct) sum Af = Ni+N-i,- the defect of the operator 
A , - has the dimension 2d = 2. Under this assumption the deficiency index of the restricted 
operator A is equal to (d, d). The defect contains elements from the domain D of A, so the 
restricted operator A = A\ Dq is not densely defined, and hence the adjoint operator does 
not exist. However, under the above condition one may use the formal adjoint Aq defined 
on defect: 

d 

u = y^[x;,j e M + xi-i e { _i] G TV, 
i=i 

by the von-Neumann formula, see for instance pq . 

d 

1=1 

We will use the formal adjoint during the construction of extensions (see 0, E3|) since, 
for closed operators with the finite deficiency indices, the construction of the extension is 
actually developed in the defect. In order to use the symplectic version of the operator- 
extension techniques we introduce in the defect M a new basis w± t t, on which the formal 
adjoint Aq is correctly defined due to the above assumption of non-overlapping: 

fl,i f _ A f f l,i f _ 7 « 
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AqW,, + = w lt _, A+w h _ = -w l+ . 



The following formula of integration by parts with abstract operators was proved in 
230): 



Lemma V.l Consider elements from the domain of the (formal) adjoint operator Aq : 

A 

A pi 1 pj 
V = Vo+ A^I^'A^I^- 

with symplectic coordinates 

i i 

Then, the boundary form of the formal adjoint operator is equal to 

J A (u,v) = (A+u,v) - (u,A+v) = (t + ,e_) N - (£!,fX (4°) 

The fundamental Krein formula 0, 0| for Generalized Resolvents of Symmetric Operators 
can be easily derived from the next statement proved in[z[J]. It helps to solve the non- 
homogeneous equations with the adjoint operator (A* — XI)u = f. In our situation similar 
to it is used in course of calculation of the scattering matrix , see below 



Lemma V.2 The symplectic coordinates £ G iV of the vector-valued function 

A + il -u A 1 -u 

U = A^XI ^ := Uo + A=il*+ ~ A=U*-> (41) 

which satisfies the adjoint equation [Aq — XI]u = 0, are connected by the formula 
The operator-function 

with a positive imaginary part in the upper half-plane 3 A > coincides actually with U Q- 
function", introduced by M. Krein (|54j). It serves an abstract analog of the celebrated 
function used by H. Weyl in course of construction of self-adjoint extensions of the second 
order differential operator of Schrodinger type, see, for instance j^, where the original 
construction by H.Weyl is presented. The operator-function M. exists almost everywhere 
on real axis A and has a finite number of simple poles sitting at the eigenvalues k i of A. 
This function plays an important role in description of spectral properties of self-adjoint 
extensions of symmetric operators, see 

BBS. 

The boundary form J(U, V) of the orthogonal sum Aq = Aq © Iq of the restricted 
matrix and the differential operator on elements (U, u) = U on the orthogonal sum of the 
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corresponding spaces is calculated as a sum of the forms ()39I40|) . The self-adjoint extensions 
of the corresponding operator A = A Q)lo are obtained as restrictions of the adjoint operator 
Aq~ = Aq © Iq onto Lagrangian planes of the form J(U, V). These planes may be defined 
by the boundary conditions connecting the symplectic coordinates U'(0), U(0), of 
components of the corresponding elements in the deficiency subspaces. For instance, one 
may select a d- dimensional operator (3 : E — > £ and define the Lagrangian plane Lg by the 
boundary condition 

/ TT'((W \ I R R \ I TT(ft\ \ 

(43) 

This condition defines a self-adjoint extension Ap of the restricted operators Aq © Iq in 
L 2 (R, E + ) © E A . The absolutely continuous of the operator Ap coincides with the spectrum 
of the exterior part of the model, hence - with the spectrum of the trivial component l A 
of the split operator C A . The corresponding eigenf unctions of on the minimal spectral 
band A A A A containing the Fermi-level can be found via substitution into above boundary 

max mtn J u 

condition the column combined of Anzatzes (|22|) and ()41|) 

m = 





e -iK+x u + e iK+x Su 
A+il , 



CU 

A-XI^+ 

It gives the linear equation for Scattering matrix: 

-iK + (u-Su)\ = ( f3 00 (3 01 \ ( v + Sv 

~C ) ) { Ml 

Solving this equation we can obtain the Scattered waves and the formula we needed: 

Lemma V.3 The Scattering matrix for the constructed model is an analytic function of the 
spectral parameter \: 

1 > iK + - [p 00 - (3 01 M(3 10 ] ' 1 ] 

with denominator of the fraction preceding the numerator. 

Now we will choose the parameter (3 such that the operator-function [j3 00 — /3 01 J\A/3 10 \ acting 
in E + coincides with the essential component T>J\f of the DN-map of the intermediate op- 
erator. Really, denote by Q s the spectral projection onto the eigen-space of A corresponding 

2 

to the eigenvalue k s framed by projections onto the deficiency subspace N t = N 

Q s =Re s )(Re s . 
Then the above expression takes the form : 

[(3 00 - f3 01 Mf3 10 ] = 



Ep^WA. (45) 
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We will define the boundary parameters /3 10 , /3 01 = j3^ Q later, but once they are de- 
fined, we choose (3 00 such that the first summand in the right side of (f43|) vanishes: 
Aw + Yl k s f3 01 Q g f3 10 = 0. Then the Scattering matrix take the form: 



iK - V ' 

m = + ^q-. , (46) 

, which coincides with the essential Scattering matrix if and only if the corresponding Krein 
-function 

E^rj^QAo (47) 



8=1 



coincides with with the essential part XW T of the DiV-map of the intermediate operator 
the finite sum of polar terms 



T „ dtp a \ / p 9(^5 



Summarizing above results we obtain the following conditional statement statement for the 
extension constructed based on the boundary condition (|43|) in case when N i D iV _. = or 
dim E . = 1: 



Theorem V.l TTie constructed operator A g zs a solvable model of the Quantum network 
on the essential interval of energy if and only if the dimension of the space E A coincides 

2 

with the number of eigenvalues on the essential interval of energy, the eigenvalues k of the 

N T 2 

operator A = ^ _ ^ s e „) ( e s > with some ortho-normal basis e s , coincide with eigenvalues of 
the intermediate operator on essential interval of energy and there exist a subspace N. =N 

and operator /3 01 : N i — > E such that for that ortho-normal basis {e s } _ in E A 



P +^t = ^P.e a (49) 



for all s, s = 1, 2, 3, . . . , N T . 



In case when only one resonance eigen-value k Q of the intermediate operator sits on the 
essential spectral band the obtained model Scattering matrix 



S(P) = ^ (50) 

iK + + ^(3 01 Q (3 W 

is a one-pole approximation of the Scattering matrix of the network. The condition of 
the above theorem are obviously fulfilled for one-pole approximation, when P + -^ ^ 0, 
d — 1, N T = 1, and (3 is one-dimensional operator mapping the one-dimensional subspace 
N i onto the resonance entrance subspace in E + spanned by P + -p L . For thin or shrinking 



34 



networks one can estimate the deviation of the one-pole approximation from the exact 
Scattering matrix on the network based on material of section 3. It will be done elsewhere. 

Note that there exist an intrinsic connection between the matrix model of the network we 
propose and method of derivation of the asymptotic of different series of eigenf unctions of 
Laplace operator on a composite domain with shrinking channels : both are actually based 
on a shape of eigenfunctions in limiting domain. We hope to develop a comparative analysis 
of both problems in a forthcoming publication. 

The above one-pole approximation has also another important role in perturbation theory: 
as a sort of a "jump-start" in analytic perturbation procedure. This role may be clarified by 
the following observation. Consider a simplest network with only one semi-infinite wire and 
the Fermi level sitting on the first spectral band multiplicity 1. Then the component of the 
corresponding solvable model in the open channel is presented by the Schrodinger equation 

—u" = Xu := p u, < x < oo, 

with one-dimensional entrance subspace E + . Construct the solvable model of the network 
assuming that the "inner Hamiltonian" A in course of construction of the solvable model 
is restricted to a symmetric operator with one-dimensional deficiency subspace N. := N 
spanned by the vector e, and /3 01 := (3 is the map of N onto E + . Then the model Scattering 
matrix is scalar : 

iP ~ (I* E s ^AQs 
S„(P) = y— (51) 

ip + P E fl 

with Q s = | (e, e s ) | 2 . Zeroes k s ((3) of the Scattering matrix (JoTj) - resonances - sit in upper 
half-plane and approach the points =tfc s , k_ a = —k B , when f3 — > 0. They can be found 
approximately in form of branching continuous fractions, see (96| as solutions of a chain of 
algebraic equations. For instance the resonance p = k (/3) created from the eigenvalue k Q of 
the operator A ( more precise : from the point +k ) can be obtained as a solution p = k (/3) 
of the equation 

P = K l f{1 + K>Q °„. T . (52) 



2 

Another resonance created at the point — k from the eigenvalue k sits at the symmetric 
point —k {(3) with respect to the imaginary axis. Other resonances k a (/3) can be found from 
similar equations with ±k substituted by ±fc s . All functions k s (f3) are analytic functions 
of P, k in small neighborhoods of (0, ±& s ). They sit in the upper half-plane symmetrically 
with respect to the imaginary axis k_ a (f3) = —k s (/3). The Scattering matrix (jHTj) is unitary 
on the real axis k and has poles at the complex-conjugate points k a (0) in the lower half-plane 
, hence it is presented by the finite Blashke product which tends to 1 when \k\ —>■ oo: 

*«-n;-5ig- 

The outer component of the scattered wave is presented as 

^ = e~ ipx + S\k)e- ipx , 
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and fulfills proper boundary condition at the place of contact with the model quantum dot. 

Though the resonances depend analytically of (f3, p), neither Scattering matrix (}5Tj) 
nor the Scattered wave depend analytically of (/3,p) on the product (Up,C p ) of a small 
neighborhood of the origin in (3 - plane and the whole p-plane. The analyticity is lost due 
to presence of points ±k a where the resonances are created at the origin of the /?-plane. But 
it is present for "separated" factors of the Scattering Matrix. Consider for small (3 two non 
intersecting subsets K and K° of the set K = {k s } and the corresponding factors S^, S° of 

the Scattering matrix with resonances created in K , K° respectively: 

P~ k s oV n TT P~K 



s » = n Et- s » = n 



p — K a H p 



Here the first factor is an analytic function of p) when (3 is small and p sits near K° , and 
, vice versa, the second factor is analytic for small f3 and p sits in a small neighborhood of 
K . In particular the first factor can be expand into Taylor series in ((3, p) for small (3 G Up 
if p sits in some small neighborhood U° of K° , and the second factor is presented by the 
Taylor series in ((3, p) if (3 £ Up and p sits in a small neighborhood U of K . On the other 
hand, the first factor is not an analytic function of (f3, p) if (3 e U„ and p is near to K . 
This observation may be interpreted as a fact of the perturbation theory the following way. 

Assume that k > and consider the representation of the Scattering matrix (J53j) in form 
of product of two factors : the product containing the resonances k ((3), — k (/3) created 
from the eigenvalue k Q 

s f>^ _ \p-KiP)\ \p + K 



[p-k (P)][p + k ((3)Y 
and the complementary factor 

One can derive from the above dispersion equation ()52|) that the resonances are analytic 
functions of the perturbation parameter (3 near the origin f3 = 0, but the factor (p) is not 
for the reason discussed above. 

Theorem V.2 There exist a one- dimensional perturbation C 1 ^ of the operator 



Lu = —u". u 



"o 







with a non-trivial inner component, such that the scattering matrix of the pair (C , l ) 
coincides with — S^(p). Then the Scattering matrix S of the complementary pair {C , C 13 ) 
is equal to the complementary factor —S° (p): 



S» = S» S>). 

The complementary factor is an analytic function of ((3, k) on the product of a small neigh- 
borhood of the origin in (3-plane and a small neighborhood K of (k , —k ) in p-plane. 
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2 

Proof Consider the one-dimensional operator A with a positive eigenvalue k . We associate 

2 

with it the Krein function M. = A . Then restrict the operator l onto smooth functions 

vanishing at the origin x — 0, and construct the self-adjoint extension A Q defined by the 
hermitian matrix {(3} which connects the boundary values £ + , — M.^ + in i?^ with boundary 
values tf(0), tf'(0) in E + : 



(54) 



Substituting into the above equation the corresponding Ansatz for the wave-function in the 
outer space 





we obtain an explicit expression for the corresponding Scattering matrix S in terms of the 

4 

inner Krein-function M.{\) = and the boundary parameters: 



S(A) 



ip + 




\P 01 \ 2 M 
1+P U M 


ip — 


Poo 


\I3 01 \ 2 M] 



Select the parameter f3 u such that 1 — /5 n /c =0, then 1 + (3 1X A4 = /3 n z K x , and the 

2 

Scattering matrix takes the form with A = p 



S(A) 



Comparing this expression with 



ip + p 



\P i I 



-j 1 2 2 

Pm \ K 2 
— p 



/3 n (l+« ) (l+« 



1^01 I' 



I^OjjV 2 



-S(P) 



[(p-*o(0)][P + *o(0)] 



[p-fc o ( / 3)][p + ^( / 3)] 

V+2^A; (/j)p+|A; (^)| 2 
p 2 -2^A; (/%-|W)| 2 1 h 
we obtain the equations for selection of parameters of the extension : 



Am 



I0oi I 



0ii (1+k ) 



1 



i/5 i r« a 



|A; (/?)| 2 2%k Q ((3) (l + « 4 )/3n 

Solution is not unique, but may be parametrized by the functional parameter /5 , 

1 



ft 



2/3 M 3fc (/3) - |*oOT 



(55) 
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In particular we may choose: f3 00 = \k {f3)\ 2 (Ssk (f3)) . Then 

1 

K "W 

To finalize the calculation we notice that both \k (/3)\ 2 and 5sk (/3) are real analytic functions 
of /3, hence k 2 is an analytic function of (3. The boundary parameter /3 00 has a pole at the 
origin since k o (0) — > k (0) = k (0) when j3 — > 0. Other boundary parameters may be 

4 

also chosen based on above equations (J55J) as /3 U = \k \ 2 , \j3 01 \ 2 = ■ Alternative 

44 
choice is = |fc | 2 , /3 00 = „ , \(3J 2 = J.^°J 4 , Ai = l fc ol~ 2 ; with K ■= K(P)- This 

accomplishes the construction of the operator A^. 

The operator A^ permits to make a sort of a "jump-start" in perturbation proce- 
dure for the pair of operators A^}: the scattering matrix S (k) for the pair 

I (I © H) , A^ I is calculated in explicit form, and the complementary scattering matrix for 

the pair |a^, A^ j, which coincides with S°, is an analytic function of (ft, p) near to the 
resonance k and can be calculated via standard analytical perturbation procedure. 

Though the "jump-start" can help to develop the analytic perturbation procedure, ac- 
cording to the preceding result, but exact calculation of the pole k o {0) and corresponding 
projection for the perturbed operator A^, which were used in course of construction of the 
"jump-start", may be a tricky analytical problem, comparable with the original spectral 
problem. On the other hand, the asymptotic expansion for these data for small f3 may be 
useful for obtaining estimates, but does not help to derive analyticity of the complementary 
factor of the Scattering matrix constructed for the corresponding approximate "jump start" 
factor, constructed based on approximate asymptotic data. 



VI. CONCLUSION: INTERMEDIATE OPERATOR IN ANALYTIC PERTURBA- 
TION TECHNIQUES 

Classical analytic perturbation techniques, see for instance |46j], is developed for additive 
perturbations of operators with discrete spectrum A — > A + eB. For perturbation of the 
continuous spectrum Kato suggested another techniques based on Wave operators. Detailed 
study of the perturbation of operators with continuous spectrum in model situations, es- 
pecially for the Friedrichs model, see |29l |. was done in numerous mathematical papers and 
books, see and the references therein and also jiO, EH |4l[ EH • It was already clear for 
specialists long ago, that the classical techniques of the analytic perturbation theory does 
not work generally on continuous spectrum (actually it works below the threshold of creation 
of resonances, see ji^]). Probably I. Prigogine alone insisted on modification of the analytic 
perturbation techniques to cover the case of continuous spectrum. 

According to evidence from long-time Prigogine's collaborator Professor I. Antoniou, 
I. Prigogine attempted also to use the idea of an intermediate operator as a base for devel- 
opment of the analytic perturbation techniques for operators with continuous spectrum. He 
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tried to construct the intermediate operator A 1 as a function of the non-perturbed operator 
A with a hope that after the first step A —>■ A x the analytical perturbation procedure for 
the perturbed operator A + eV presented as A 1 + eV becomes convergent. This simplest 
suggestion appeared to be not so efficient as expected, and Prigogine abandoned this idea 

finally, see @, H, E3, H El 

Note that actually the idea of the intermediate operator in another context appeared 
in Splitting Method invented by I. Glazman, see jiE]- This method radically transformed 
in seventies the whole picture of the qualitative spectral analysis of singular differential 
operators, but was never used in context of the analytic perturbation procedure. 

We actually suggested revisiting the idea of the intermediate operator again in a special 
case of quantum networks where the candidate for the role of an intermediate operator 
is almost obvious. In our case it is not a function of the non-perturbed operators, but is 
obtained from it just by a finite- dimensional perturbation. Actually we succeeded to apply the 
splitting method to Quantum Networks because of the step-wise structure of the continuous 
spectrum in wave-guides. This minor change of the original idea of Prigogine appears to be 
useful for perturbation problems in wave-guides. In our case the idea of the intermediate 
operator is automatically combined with resonances, which was anticipated lon g ag o by 
Poincare and intensely discussed by Prigogine with references to Poincare, see jzJEJ 

The only essential, and even decisive difference of our version of the intermediate operator 
from the version which I. Prigogine attempted to use, is the natural principle of spectral 
locality we impose: we construct the intermediate operator and the scattered waves via the 
special perturbation procedure locally- in a certain spectral interval, - but not globally, as 
Prigogine planned. It is quite easy to define this interval mathematically for the Quantum 
Networks: it is the maximal interval containing the Fermi level, j^j, where the spectral 
multiplicity n is constant and the scattering matrix is a square matrix size n x n. This 
interval may be either a spectral band corresponding to some channel in one of semi-infinite 
wires, of a sum of several spectral bands (sort of a "spectral terrace") corresponding to 
different wires, with properly hybridized scattered waves. The structure of the corresponding 
scattered waves within this spectral interval is manifested by the structure of the scattering 
matrix and corresponds to the n-dimensional input-output. The structure of scattered waves 
on the neighboring spectral bands is different. 

It is worth to notice, that slightly modified procedure permits to construct an hierarchy 
of intermediate operators also for some operators with absolutely-continuous spectrum of 
constant multiplicity, in particular for Friedrichs model. 

In actual paper we discussed simplest wave-guiding systems composed of straight wave- 
guides (quantum wires). The convenience of this trivial environment appears from existence 
of explicit formulae for oscillatory modes in straight wave-guides: the system of thresholds 
and the structure of the entrance sub-spaces in this environment is trivial. Study of more 
interesting systems of periodic or asymptotically straight wave-guides with rapidly decreas- 
ing curvature requires a combination of our methods with Lippmann-Schwinger techniques 
developed in case of wave-guides with non-trivial geometry for discrete and especially for 

continuous spectrum, see 0, HBO SB 



On the other hand, the suggestion to recover resonances from the discrete spectrum of 
some intermediate operator can be easily applied to any wave-guiding problem in trivial 
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environment, with linear vertices, on a reasonable price of preliminary study of the "in- 
termediate" dynamics. In particular, it permits to suggest alternative approach to study 
of trapped modes in acoustics and hydro-dynamics, investigated in [3, [23 , as well as to 
classical theory of electro-magnetic wave-guides with local variation of curvature, [20]. An 
important advantage of our techniques is: the possibility of "fitting" of the constructed 
"solvable model" based of results of straightforward computing for the intermediate prob- 
lem. 
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